HIGH-RESOLUTION OBSERVATIONS OF MOLECULAR 
LINES TOWARD THE HOT CORE G28.20-0.04N 



Sheng-Li Qin^, Maohai Huang^, Yuefang Wu^, Rui Xue^, Sheng Chen^ 

'66 ; 

O ■ ABSTRACT 
O ■ 

(N ; 

W)| We present the results from arcsecond resolution observations of various line 

^ I transitions at 1.3 mm toward hypercompact HII region G28.20-0.04N. With the 

^ I SMA data, we have detected and mapped the transitions in the CH3CN, CO, 

CN I ^^CO, SO2, OCS, and CH3OH molecular hues as well as the radio recombination 

,— 1| line HSOa. The observations and analysis indicate a hot core associated with 
G28.20-0.04N. The outflow and possible rotation are detected in this region. 

6 ■ 
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C§ I and dynamics — ISM:molecules — radio lines: ISM — star: format ion 
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cn : 1. INTRODUCTION 

' The scenario of massive star formation remains unclear and has been observationally 

00 ' challenging because of large distances, clustered formation environments, and shorter evo- 
^ ■ lutionary timescales of massive stars. Ultracompact HII (UCHII) regions are considered 
O ■ signposts of massive star formation, but do not represent the earliest stage of massive star 
' forming process (e.g., Churchwell 2002). Hot molecular cores are defined as compact (< 0.1 
^ ■ pc), dense (> 10'' cm~'^) and warm (> 100 K) molecular cloud cores (Kurtz et al. 2000). The 
• observations from various wavelengths have suggested that hot cores are the sites for massive 
star formation, and that they represent the early phase of the evolution prior to UCHII re- 
gions (Kurtz et al. 2000; Gibb et al. 2000, 2004; Churchwell 2002). So far only twenty or so 
hot cores from high-resolution observations have been reported in the literature. The obser- 
vations of hot cores are important for understanding the evolutionary sequence and physical 
conditions of massive star formation. Hypercompact HII (HCHII) regions have smaller sizes 
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and higher densities compared with UCHII regions, and they probably correspond to a tran- 
sition phase from hot cores to UCHII regions. 

The HCHII region G28.20-0.04N at a distance of 5.7 kpc has been observed at radio and 
millimeter wavelengths (Fish et al. 2003; Sollins et al. 2005; Sewilo et al. 2004, 2008; Keto 
et al. 2008). Masers, rotation, inflow, and outflow in G28.20-0.04N have been revealed from 
observations of molecular and radio recombination lines at centimeter wavelengths (Argon et 
al. 2000; Menten 1991; Sollins et al. 2005; Sewilo et al. 2008). However, the star formation 
environment and physical conditions of the hot core in this region are still poorly understood. 

Due to the compact and dense nature of hot cores, high angular resolution obser- 
vations using molecular lines with high critical densities and excitation temperatures at 
(sub) millimeter wavelengths are crucial to uncovering the physical conditions and kinemat- 
ics in G28.24-0.04N. In this letter, we present the Submillimeter Array (SMA) [^observations 
of molecular lines at 1.3 mm toward the hot core in G28.24-0.04N. 



2. DATA 

The data are from SMA archive. Observations toward G28.20-0.04N were carried out 
with the SMA in September 2005, at 220 GHz (lower sideband) and 230 GHz (upper side- 
band) with a frequency resolution of 0.4125 MHz and an angular resolution of 1". The 
mean system temperature was 100 K. In the database, QSO 3c454.3 and Uranus were ob- 
served for bandpass and fiux-density calibrations. QSOs 1743-038 and 1911-201 were also 
observed for the antenna-based gain corrections. The calibration and imaging were done in 
Miriacfl. Multiple lines were detected in both sidebands. Using the task UVLIN in Miriad, 
the continuum level was determined by a linear fitting to the line-free channels. Then, the 
calibrated (u, v) data were separated for two output (u, v) data sets: the continuum and 
the continuum-subtracted spectral lines. In order to improve sensitivity, we smoothed the 
H30a, CO and ^'^CO lines to 2 km s~^ velocity resolution, and smoothed rest of the lines 
to 1 km s~^ velocity resolution. Self-calibration was performed to the continuum data. The 
gain solutions from the continuum were applied to the line data. The synthesized beam 
size of the continuum and line images with robust weighting was approximately 1.3" x 0.8" 
(P.A.=-89.6°). 



^The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the 
Academia Sinica Institute of Astronomy and Astrophysics and is funded by the Smithsonian Institution and 
the Academia Sinica. 

^http: / /sma-www. cfa.harvard.edu/miriadWWW 
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3. RESULTS 

The molecular lines were identified using the CDMS and JPL databases (Miiller et al. 
2005; Pickett et al. 1998), and the data by Sutton et al. (1985) & Nummehn et al. (1998). 
The molecular lines CH3CN, CO, ^'^CO, SO2, OCS, CH3OH, and the radio recombination 
fine H30q; were detected and identified from the lower and upper sidebands of the SMA 
observations. Table 1 fists the parameters of the detected molecular transitions. 

Combining the lower and upper sideband data, we made continuum image from line- 
free channels at 1.3 mm, as shown in Fig. 1(a). The rms (1 a) noise level is 0.003 Jy 
beam"-^. The observations show an unresolved core peaked at R.A.(J2000.0)=18''42™58.120^ 
(AR.A. = ±0.02"), decl.(J2000.0)=-4°13'57.40" (Adecl. = ±0.02") with a peak intensity 
of 0.67±0.02 Jy beam~^. The total flux density and deconvolved source size from Gaussian 
fit are approximately 0.89±0.03 Jy and 0.6" x 0.2"(P.A. = -10.1°), respectively. The thick 
contours in Fig. 1(a) are the integrated intensity map of the radio recombination line H30a. 
The SMA observations show that the peak position of the continuum is consistent (within 
0.03") with that of the H30q; line image in G28.24-0.04N. Most likely the continuum emission 
in this observation comes from free-free emission (see Fig. 10 of Keto et al. 2008). 

The line images were constructed from the maps of the continuum-subtracted spectral 
channels. Figure 1(b) shows the integrated intensity contours of the CH3CN (122 — II2) tran- 
sition. The major axis of the line image is approximately 1.4", corresponding to a projected 
linear size of <0.04 pc. In order to avoid the cancellation of the integrated intensities of the 
emission (positive values) and absorption (negative values) in the overlapping region along 
the hne-of-sight, the emission and absorption in the CO and ^^CO were handled separately 
for the integrated intensity maps, as done by Qin et al. (2008). Figure 1(c) presents the in- 
tegrated intensity map of the ^^CO. The absorption against the continuum and the emission 
were observed as shown in Fig. 1(c). No absorption was evident in the H30q; and CH3CN 
lines. Away from the continuum, the absorption was not observed. Figure 2 presents the 
spectra of the molecular line transitions. The spectra were extracted from the channel maps 
at the peak positions of the line images. In Fig. 2, the emission spectra have similar profiles, 
showing single emission peak with similar line center velocity and line width. Gaussian fits 
were performed on all the spectra. The peak intensity (Ip), full width at half maximum 
(Ay), and central line velocity (Vlsr) from the Gaussian fits are summarized in Table 1. 
The radial velocities range from 95 to 96 km s~^ (except for ^^CO, CO and H30q; lines). 
Compared with the SO2, OCS, CH3OH and CH3CN, the ^^CO line with lower excitation 
temperature shows both absorption and emission components. Absorption can be observed 
if the excitation temperature of gas in front of the background continuum is lower than the 
brightness temperature of the continuum. Absorption is also observed in the CO line. 
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4. DISCUSSIONS 
4.1. Hot Core 

Methyl cyanide (CH3CN) has been proved as an ideal probe to determine the kinetic 
temperature and column density of molecular gas (e.g., Remijan et al. 2004). Seven K- 
components of the CH3CN {J — 12 — 11) transition were detected by the SMA observations. 
From the integrated intensities, in the limits of optically thin and local thermodynamic 
equilibrium (LTE), the rotation temperature and column density are estimated using a ro- 
tation temperature diagram (Goldsmith & Langer 1999; Liu et al. 2002). Figure 3 shows 
the rotation temperature diagram. A linear least-square fit is performed toward the seven 
CH3CN transitions. The derived rotation temperature and beam-averaged column density 
are 308±22 K and (1.6±0.3) xlO^^ cm~^, respectively. In the Orion molecular cloud, the 
CH3CN transitions were only detected in the hot core and the compact ridge, giving frac- 
tional abundances relative to H2 of 7.8x10"^ and 3.2x10"^'^, respectively (Blake et al. 1987). 
Taking source size of 0.04 pc and the fractional abundance in the Orion cloud, we inferred the 
H2 density of > 10^ cm~'^ in G28.24-0.04N. The relatively higher gas temperature (308±22 
K), H2 density (> 10^ cm~^) and smaller size (<0.04 pc) indicate a hot core in this region. 

Following Gibb et al. (2000), we adopt 3 x 10^^ cm~^ as a typical H2 column density. 

Using this value, we calculate a fractional abundance of CH3CN relative to H2 of 5 x 10~^. 
The derived abundance is close to that in the Orion hot core (Blake et al. 1987). The 
derived rotation temperature agrees (within 2a of the least-square fit) with the vahie of 280 
K estimated by the NH3 lines (Soillns et al. 2005). The similar gas temperatures derived 
from the CH3CN and NH3 lines suggest a close relationship between the two N-containing 
species. The gas-phase chemistry likely dominates the pathways to produce CH3CN (Rodgers 
& Charnley 2001), in which NH3 evaporates from grain surfaces and CH3CN is then formed 
via gas-phase reactions at high temperatures. Based on the high temperature (300 K) 
gas-phase chemical model (Rodgers & Charnley 2001), if NH3 is injected for the chemical 
reactions, the derived fractional abundance of the CH3CN of 5x10^*^ corresponds to a time- 
scale of 1.5x10'^ yr which is in good agreement with those (1.9x10'^ — 5.7 x lO'' yr) observed 
in other hot cores (Kurtz et al. 2000). The CH3CN, SO2, OCS and CH3OH spectra have 
similar radial velocities and line widths, suggesting that the N-bearing, 0-bearing and S- 
bearing molecules probably originated from the hot core. 
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4.2. Kinematics 

For the ^'^CO absorption component in Fig. 1(c), the derived optical depth of the hne 
peak is 0.57 by use of the hne to continuum ratio (cf. equation (1) of Qin et al. 2008). The 
intensity of 1 Jy beam^^ in this observation corresponds to a brightness temperature of 26 
K under the Rayleigh- Jeans approximation. Given the beam fiUing factors of the hne and 
continuum of 0.5, the upper hmit of the excitation temperature in the ^^CO hne is 35 K (cf. 
equation (5) of Qin et aL 2008). Assuming the gas is in LTE and the emission component 
(with integrated flux density of 9.1 Jy beam~^ km s~^) is optically thin with abundances 
of [CO]/[^3co]=89 and [CO]/[H2]=lxlO-^ the H2 column densities of the absorption and 
emission components are ~ 4.6x10^^ and 6.8x10^^ cm~^, respectively. The major axes of the 
absorption and emission components are approximately 0.025 and 0.04 pc. The estimated H2 
densities and masses are 5.9x10^ cm~^, 1.5 Mq and 5.5x10^ cm~'^, 1.9 Mq for the absorption 
and emission components, respectively. Compared with the rotation temperature and H2 
density traced by CH3CN hnes, the lower excitation temperature and H2 densities suggest 
that the ^^CO is located outside of the hot molecular core. The average Vlsr from multiple 
lines is an approximation of the systemic velocity of the molecular cloud (Sutton et al. 1991). 
The average Vlsr of the emission lines (except for H30a, CO and ^^CO) is 95.4±0.5 km s^^ 
(see Table 1). We adopt 95.4±0.5 km s~^ as the systemic velocity of the cloud. The spectrum 
in the ^"^CO line (see Fig. 2(j)) shows that the absorption is blue-shifted and the emission 
is red-shifted with respect to the systemic velocity, consisting of a P Cygni-like profile. A 
NE-SW velocity gradient and broad line wings in NW-SE direction have been observed in 
the radio recombination lines (Sewilo et al. 2008). The authors proposed a rotating motion 
responsible for the observed velocity gradient and an outflow along the rotational axis in 
the NW-SE direction. In Fig. 1(c), the peaks of the absorption and emission are ahgned 
with the rotational axis. The spectral profile along with the ^^CO distribution appear to 
indicate an outflow from the molecular core, where the blue-shifted absorption (in front of 
the continuum) is flowing toward us and the red-shifted emission (behind the continuum) is 
moving away in the opposite direction. The CO spectrum shows a complex with red-shifted 
emission and both blue and red-shifted absorption with respect to the systemic velocity (see 
Fig. 2(i)). The infall component in this region were observed by red-shifted absorption in 
the NH3 line (Sollins et al. 2005). Likely the blue-shifted outflow traced by the CO hne in 
this observation is mixed with the infall. Higher spatial and spectral resolution observations 
are needed to verify the kinematics traced by the CO hne. 

The position-velocity diagram across the peak of the continuum along the NE-SW di- 
rection is constructed from the OCS hne (see Fig. 4). In Fig. 4, the velocities of the two 
emission peaks are 95 and 97 km s~^ with 0.2" separation, respectively, indicating a velocity 
gradient in NE-SW direction. Assuming a rotating motion along NE-SW with a rotation 
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axis in the NW-SE direction (Sewilo et al. 2008) and an equilibrium between rotational 
and gravitational forces, the dynamical mass responsible for the rotation can be estimated 
by M = X r/G, where V and r are the velocity difference and spatial separation of the 
emission peaks, respectively; G is the gravitational constant. The derived dynamical mass 
is 25 Mq, which is consistent with the mass derived by Sewilo et al. (2008). 

We thank the SMA staff for making the observations possible. We thank the anonymous 
referee for constructive comments on the paper. Y.F. Wu acknowledges the support from 
NSFC under grant 10733030. 
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Table 1. Molecular Line Parameters 



Molecule 


Transition 


Frequency 


Eu 


Ip 


AV 


Vlsr 


Channel rms 






(MHz) 


(K) 


(Jy beam~^) 


(km s"-*^) 


(km s~^) 


(Jy beam^^) 


CO" 


2-1 


230538.00 


17 


-0.2±0.07 


9.9±3.0 


79.9±1.3 


0.07 










-0.3±0.07 


10.3±4.5 


94.1±1.6 












0.28±0.09 


7.1±2.1 


103.1±1.3 




13qO» 


2-1 


220398.68 


16 


-0.29±0.05 


4.8±0.8 


77.1±0.4 


0.05 










0.3±0.05 


6.3±0.9 


104.3±0.4 




OCS 


19-18 


231060.98 


111 


1.1±0.06 


3.7±0.2 


95.4±0.1 


0.06 


SO2 


lli,ii-10(),io 


221965.21 


60 


1.5±0.07 


4.1±0.2 


95.7±0.1 


0.07 


CH3OH A 


IO2-93 


231281.10 


166 


0.3±0.07 


4.6±0.8 


95.5±0.3 


0.07 


CHaCN 


I20-II0'' 


220747.26 


69 


l.liO.08 


3.8±0.3 


95.1±0.1 


0.08 




I21-II1'' 


220743.01 


76 


1.2it0.08 


4.9±0.4 


95.0±0.1 


0.08 




I22-II2 


220730.26 


97 


0.9±0.08 


4.1±0.3 


95.3±0.1 


0.08 




I23-II3 


220709.02 


133 


1.0±0.08 


4.2±0.3 


95.4±0.1 


0.08 




I24-II4 


220679.29 


183 


0.6±0.08 


4.8±0.5 


95.4±0.2 


0.08 




I25-II5 


220641.09 


247 


0.7±0.07 


3.9±0.3 


95.3±0.1 


0.07 




12fi llfi 


220594.43 


325 


0.6±0.08 


3.3±0.4 


95.8±0.2 


0.08 


ll:-!0n' 




2:-!l<)0().!)(i 




l).!)±0.2 


2I).I)±0.() 


!)2.5±0.2 


O.IKi 



"Negative values of the intensity {Ip) in the CO and ^"^CO spectra indicate the absorption. Positive values indicate 
the emission. 

''The two CH3CN transitions (12o-llo) and (12i-lli) were blended. 
°H30a line parameters are cited from the paper by Keto et al. 2008. 
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Fig. 1. — Continuum and spectral images, (a): The continuum (thin contours) is superim- 
posed on the H30q; radio recombination hue image (thick contours). The synthesized beam 
is 1.3" X 0.8", P.A.=— 89.6° (lower-right corner). The rms (1 a) noise level of the continuum 
is 0.003 Jy beam-^ The contours are -4, 4, 8, 16, 32, 64, 128, 190, and 220 a. The cross 
symbol indicates the peak position of the continuum source. The contours of the HSOo; inte- 
grated intensity are -1.24, 1.24, 3.68, 6.13, 8.58, 11.03, 13.48, 15.93, 18.38, 20.83, and 23.28 
Jy beam~^ km s~^; (b): The CH3CN (122 — II2) integrated intensity contours. The levels are 
-0.74, 0.74, 1.23, 1.72, 2.21, 2.70, 3.19, 3.68, 4.17, and 4.66 Jy beam"^ km 8"^ (c): The ^^CO 
(2-1) integrated intensity contours. The levels are —2, —1.8, —1.6, —1.4, —1.2, —1, —0.8, —0.6, 
1, 1.3, 1.6, 1.9, 2.2, 2.5, 2.8, 3.1 Jy beam~^ km s~^; the positive and negative values indicate 
the emission and absorption, respectively. 
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(o) CH3CN K=0.1 




85 90 95 100 




(I) CO 




^^^^^ 






Velocity (km s ^) 

Fig. 2. — Molecular spectra averaged over one beam. The solid and dashed curves are the 
observed spectra and the Gaussian fitting to the spectra, respectively, (a), (b), (c), (d), (e) 
and (f) are for the various K components in the CH3CN (12—11); (g), (h), (i) and (j) are for 
OCS (19-18), CH3OH (IO2 -93), CO (2-1) and ^^CO (2-1) transitions, respectively. The 
spectra are Manning smoothed for better signal to noise ratios. Negative and positive values 
of the intensity in the CO and ^^CO spectrum indicate absorption and emission, respectively. 
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Fig. 3. — Rotation temperature diagram for the observed CH3CN transitions. The hnear 
least-square fit (indicated by sohd hne) gives a rotation temperature of 308±22K. The vertical 
bars indicate the errors of In (Nu/gu) transferred from the integrated intensities. 
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Fig. 4. — Position-velocity diagram cutting along NE-SW direction across the continuum 
peak. The diagram is constructed from the OCS transition. The contours are 0.19, 0.38, 
0.56, 0.75, 0.94 and 1.13 Jy beam-^ 



